Search for: All records

We present a comprehensive quantum mechanical study of stereodynamic control of HD + He and D2 + He collisions that have been probed experimentally by Perreault et al. [J. Phys. Chem. Lett. 13, 10912 (2022)] using Stark-induced adiabatic Raman passage (SARP) techniques. Our calculations utilize a highly accurate full-dimensional H2 + He interaction potential with diagonal Born–Oppenheimer correction appropriate for HD and D2 isotopomers. The results show that rotational quenching of HD from j = 2 → j′ = 0 in v = 2, j = 2 → j′ = 1 in v = 2 and v = 4, and j = 4 → j′ = 3 in v = 4 is dominated by an l = 1 shape resonance located between 0.1 and 1.0 cm−1. For collision energies less than 0.1 cm−1, isotropic scattering prevails. An l = 1 resonance centered around 0.02 cm−1 is also found to dominate the j = 2 → j′ = 0 and j = 4 → j′ = 2 transitions in v = 4 for He–D2 collisions consistent with our prior studies of Δj = −2 transition in He + D2(v = 2, j = 2) collisions. Our analysis does not support the hypothesis of Perreault et al. [J. Phys. Chem. Lett. 13, 10912 (2022)] that a strong l = 2 resonance controls the angular distribution for Δj = −2 transition for both systems. Despite improvements in the development of the potential energy surface, a good agreement with SARP experiments for v = 2 is achieved only when contributions from collision energies less than 1.0 cm−1 were excluded in the computation of velocity averaged differential rate coefficients for both systems. This could be due to some uncertainties in the velocity spread in the experiment that employs co-propagation of the collision partners and possibly, the neglect of transverse velocities in the simulation of the experiment. 

Symmetry-adapted perturbation theory (SAPT) is a popular and versatile tool to compute and decompose noncovalent interaction energies between molecules. The intramolecular SAPT (ISAPT) variant provides a similar energy decomposition between two nonbonded fragments of the same molecule, covalently connected by a third fragment. In this work, we explore an alternative approach where the noncovalent interaction is singled out by a range separation of the Coulomb potential. We investigate two common splittings of the 1/r potential into long-range and short-range parts based on the Gaussian and error functions, and approximate either the entire intermolecular/interfragment interaction or only its attractive terms by the long-range contribution. These range separation schemes are tested for a number of intermolecular and intramolecular complexes. We find that the energy corrections from range-separated SAPT or ISAPT are in reasonable agreement with complete SAPT/ISAPT data. This result should be contrasted with the inability of the long-range multipole expansion to describe crucial short-range charge penetration and exchange effects; it shows that the long-range interaction potential does not just recover the asymptotic interaction energy but also provides a useful account of short-range terms. The best consistency is attained for the error-function separation applied to all interaction terms, both attractive and repulsive. This study is the first step toward a fragmentation-free decomposition of intramolecular nonbonded energy. 

An impressive collection of accurate two-body interaction energies for small complexes has been assembled into benchmark databases and used to improve the performance of multiple density functional, semiempirical, and machine learning methods. Similar benchmark data on nonadditive three-body energies in molecular trimers are comparatively scarce, and the existing ones are practically limited to homotrimers. In this work, we present a benchmark dataset of 20 equilibrium noncovalent interaction energies for a small but diverse selection of 10 heteromolecular trimers. The new 3BHET dataset presents complexes that combine different interactions including π −π, anion−π, cation−π, and various motifs of hydrogen and halogen bonding in each trimer. A detailed symmetry-adapted perturbation theory (SAPT)-based energy decomposition of the two- and three-body interaction energies shows that 3BHET consists of electrostatics- and dispersion-dominated complexes. The nonadditive three-body contribution is dominated by induction, but its influence on the overall bonding type in the complex (as exemplified by its position on the ternary diagram) is quite small. We also tested the extended SAPT (XSAPT) approach which is capable of including some nonadditive interactions in clusters of any size. The resulting three-body dispersion term (obtained from the many-body dispersion formalism) is mostly in good agreement with the supermolecular CCSD(T)−MP2 values and the nonadditive induction term is similar to the three-body SAPT(DFT) data, but the overall three-body XSAPT energies are not very accurate as they are missing the first-order exchange terms.